Realization of Arbitrary Gates in Holonomic Quantum Computation 
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Among the many proposals for the reahzation of a quantum computer, holonomic quantum com- 
putation (HQC) is distinguished from the rest in that it is geometrical in nature and thus expected 
to be robust against decoherence. Here we analyze the realization of various quantum gates by 
solving the inverse problem: Given a unitary matrix, we develop a formalism by which we find loops 
in the parameter space generating this matrix as a holonomy. We demonstrate for the first time 
that such a one-qubit gate as the Hadamard gate and such two-qubit gates as the CNOT gate, the 
SWAP gate and the discrete Fourier transformation can be obtained with a single loop. 
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I. INTRODUCTION 

Quantum computing is an emerging scientific disci- 
pline, in which the merging and mutual cross-fertilization 
of two of the most important developments in physical 
science and information technology of the past century — 
quantum mechanics and computing — has resulted in an 
extraordinarily rapid rate of progress of interdisciplinary 
nature. Interesting problems to address in this context 
include fundamental questions as to what are the ulti- 
mate physical limits of computation and communication. 
For introductions to quantum computing and quantum 
information processing see, e.g., Refs. |[ |j. 

Holonomic quantum computation (HQC) was first sug- 
gested by Zanardi and Rasetti in Ref. Ml. The con- 
cept has been further developed in Refs. 0, |[ 
The suggestion is very intriguing itself; quantum-logical 
operations are achieved by driving a degenerate system 
around adiabatic loops in the parameter manifold. The 
resulting gates are a generalization of the celebrated 
Berry phase to encompass a degenerate system. 
These are, in fact, non-Abelian holonomies. Due to the 
geometric nature of these gates, quantum information 
processing is expected to be fault tolerant. For instance 
the issue of timing and the lack of spontaneous decay 
are definite strengths of HQC. Here we study the con- 
struction of holonomic quantum-logic gates numerically 
for the first time via solving a certain inverse problem. 
Namely, we find the loop 7 corresponding to the desired 
unitary operator U by solving a high-dimensional opti- 
mization task. 

The paper is organized as follows: In Section we 
present the physical and mathematical background un- 
derlying our approach. Sections [II, IV and ^ comprise 
the main part of the present paper. Loop parameter- 
izations for one- and two-qubit gates are presented in 
Section HI. The numerical method is introduced in Sec- 



tion rV. Then the optimal realization of a unitary gate 
as a holonomy associated with a loop in the parameter 
space is investigated numerically in Section 0. Section VI 
discusses the results. 



II. HAMILTONIAN AND HOLONOMY 

Here we first review the concept of non-Abelian holon- 
omy to establish notation conventions. Let us consider 
a family of Hamiltonians {Hx}. The point A, continu- 
ously parameterizing the Hamiltonian, is an element of a 
manifold A4 called the control manifold and the local co- 
ordinate of A is denoted by A' (1 < i < m = dimA^). It is 
assumed that there exists only a finite number of eigen- 
values £i;(A) (1 < A: < i?) for an arbitrary X G A4 and 
that no level crossings occur. Suppose the n-th eigen- 
value e„(A) is g„-fold degenerate for any A G and 
9n = N ■ The degenerate subspace at A is denoted 
by 7i„(A). Accordingly, the Hamiltonian is expressed 
as an A^ X A^ matrix. The orthonormal basis vectors of 
7i„(A) arc denoted by {Ina; A)}; 



Hx\na; A) = £„(A)|7ia; A), {na; \\m(3; A) = 5„inSap- 

Note that there are U{gn) degrees of freedom in the 
choice of the basis vectors {\na; A)}. 

Let us now assume that the parameter A is changed 
adiabatically. We will be concerned with a particular 
subspace, say the ground state Hi (A) and we drop the 
index n to simplify the notation. Suppose the initial state 
at i = is an eigenstate |V'a(0)) = \a; A(0)) with the en- 
ergy e = possibly through shifting the zero-point of 
the energy. In fact, we are not interested in the dynami- 
cal phase at all and hence assume that the eigenvalue in 
this subspace vanishes for any \ E M. The Schrodinger 
equation is 
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whose solution may be assumed to take the form 
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iv-aW) = El/3;^W)^^/3"W• 
/3=l 



(2) 



The unitarity of the matrix Upait) foUows from the nor- 
mahzation of lipait))- By substituting Eq. into (|^), 
one finds that Uf^a satisfies 
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7;A(0)C/^a. (3) 



The formal solution may be expressed as 
U{t) = Texp AiT)dT 



AiT)dT 



f dT f dT'A{T)A{T') 

Jo Jo 



(4) 



where T is the time-ordering operator and 



Ap^{t) = {f3;X{t) 



a; X{t) 



Let us introduce the Lie-algebra-valued connection 

^ a;X{t) 



= ( /3; X{t) 



through which U{t) is expressed as 



U{t) ^ V exp 



AidX' , 



(5) 



(6) 



A(0) 



where V is the path-ordering operator. Note that Ai is 
anti-Hcrmitian, ^| — —Ai- 

Suppose the path X{t) is a loop 7(t) in Ai such that 
7(0) = 7{T) = Xq. Then it is found after traversing 7 
that one ends up with the state 



\MT)) - IMm^o^iT) (7) 

/3=1 

where use has been made of the definition lipplO)) — 
1/3; Aq). The unitary matrix 

EE U{T) = V exp ^- (j) A^dY^ 



(8) 



is called the holonomy associated with the loop 7(i). 
Note that U-y is independent of the parameterization of 
the path but only depends upon its geometric image in 
M. 

The space of all the loops based at Aq is denoted 
Lx„iM) - {7 : [0,T] ^ MhiO) = 7(T) = Aq}. (9) 



The set of the holonomy 



(10) 



has a group structure and is called the holonomy 
group. It is clear that Hol(^) C U{g). The connection 
A is called irreducible when Hol(^) = U{g). 



III. THREE-STATE MODEL AND 
QUANTUM-GATE CONSTRUCTION 

A. One-qubit gates 

To make things tractable, we employ a simple model 
Hamiltonian called the three-state model as the basic 
building block for our strategy. This is a 3-dimensional 
Hamiltonian with the matrix form 





i/Ao =e|2)(2|= (000 




(11) 



The first column (row) of the matrix refers to the auxil- 
iary state 1 2) with the energy e > while the second and 
the third columns (rows) refer to the vectors |0) and |1), 
respectively, with vanishing energy. The qubit consists 
of the last two vectors. 

The control manifold of the Hamiltonian ( pT|) is the 
complex projective space CP^. This is seen most di- 
rectly as follows: The most general form of the iso- 
spectral deformation of the Hamiltonian is of the form 
= W^HxgW^, where G U{3). Note, however, 
that not all the elements of U{3) are independent. It is 
clear that Hj is independent of the overall phase of , 
which reduces the number of degrees of freedom from 
U{3) to U{3)/U{1) = SU{3). Moreover, any element of 
SU{3) may be decomposed into a product of three SU{2) 
matrices as follows 



W-y = 



/3i ai 0> 
~ai Pi 
1, 



h 012' 
1 

-"2 (32. 



U2 



'10 

fis "3 

^0 -ag (3^^ 

U3 



(12) 

which is know as the Givens decomposition. Here aj — 



e«0j sin 6j and Pj — e^^^ cos ( 



'j . It is clear that H-y is inde- 



pendent of U3 since [Hxo,!/^] = 0. This further reduces 
the physical degrees of freedom to SU{3)/SU{2) = S^. 
The product U1U2 contains six parameters while is 
five-dimensional; there must be one redundant parameter 
in UiU2- This parameter is easily found out by writing 
the product explicitly. The result depends only on the 
combination 02 ~ ip2 and not on individual parameters. 
Accordingly, we may redefine 02 as 02 — ip2 to eliminate 
"02 • Furthermore, after this redefinition we find that the 
Hamiltonian depends only on 0i — ipi and 02 — "01 and 
hence ■0i may also be subsumed by redefining 0i and 02 , 



3 



which reduces the independent degrees of freedom down 

to cp2 ^ sys\ 

Let [z^, z'^, z^] be the homogeneous coordinate of CP^ 
and (1, ^1, ^2) be the corresponding inhomogeneous coor- 
dinate, where ^1 = z^jz^^^ = z^jz^ in the coordinate 

neighborhood with ^ 0. If we write = r^e*'^'', the /g q\ 

above correspondence, i.e. the embedding of CP^ into = ( q q ] > (1^) 

f (3), is expHcitly given by 9^ = tan"-^ rfe and ^fe = ifk- ^ ' 

The connection coefficients are easily calculated in the 
present model and are given by 

- Vsin^2e^(^=-^^) J ' ^^'^> 
I 



—I sm 



.—^ sin 2^1 sin 626 



-|sin26ii sin6'2eH'*'i 
;(<^2-<Ai) i sin=^ 6*2 sin^ 6I1 



(15) 







-isin^6l2 



2« > 



(16) 



where the first column (row) refers to |0) while the second 
one refers to |1). Using these connection coefficients, it is 
possible to evaluate the holonomy associated with a loop 
7 as 



(17) 



V exp ^- j) {Ae^dOi + Ae2d62 + A^^dcpi + ^02*^02)^ • 

Now our task is to find a loop that yields a given unitary 
matrix as its holonomy. 



B. Two-qubit gates 



Let us consider a two-qubit reference Hamiltonian 



H 



2-qubit 



73 73 o 77; 



An' 



(18) 



where 77"'*' are three-state Hamiltonians and 73 is the 
3x3 unit matrix. Generalization to an arbitrary N- 
qubit system is obvious. The Hamiltonian scales as 3''^, 
instead of the 2^ in the present model. It is also pos- 
sible to consider a model with g-dcgcnerate cigcnstatcs 
with one auxiliary state having finite energy. This model, 
however, has a difficulty in realizing an entangled state, 
without which the full computational power of a quan- 
tum computer is impossible. 

We want to maintain the multipartite structure of the 
system in constructing the holonomy. For this purpose. 



we separate the unitary transformation into a product 
of single-qubit transformations (W^ (E) W^) and a purely 
two-qubit rotation which cannot be reduced into 

a tensor product of single-qubit transformations. There- 
fore, we write the iso-spectral deformation for a given 
loop 7 as 



jy2-qubit 



^2-qubit ^ ;y6)Jj2-q"bit 
X(W^»0 W^)tW^2-qubitt_ 



(19) 



The advantage of expressing the unitary matrix in this 
form is easily verified when we write down the connection 
coefficients for the one-qubit coordinates. Namely, the 
two-qubit transformation does not affect the one-qubit 
transformation at all; 



= (a;\ 



d 

W-v 



/3;A), 



where 7' denotes a one-qubit coordinate. 

There is a large number of possible choices for VK^"''"'^'*, 
depending on the physical realization of the present sce- 
nario. To keep our analysis as concrete as possible, we 
have made the simplest choice 

^2-qubit ^^^=gi£|ll>(ll| (20) 

for our two-qubit unitary rotation. Let 
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be a two-qubit Hamiltonian before is applied. Then after the application of IVj to H'^ we have the full Hamiltonian 



^2-qubit 



ft'' 

"21 




'''3 


hi2 








hl3 








ftll+/l22 


''23 
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''33"^ ''11 




ft53e-'« 
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''32 
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/l^2e'« 


^33+^^33 



(21) 



r 



As for the connection, we find 













































icos^ 


9? 



(22) 



where the columns and rows are ordered with respect 
to the basis {|00), |01), |10), |11)}. It should be apparent 
from the above analysis that we can construct an arbi- 
trary controlled phase-shift gate with the help of a loop 
in the (6*^, ^)- or (6*^, $)-space. Accordingly, this gives the 
CNOT gate with one-qubit operations, as shown below. 



C. Some Examples 

Before we proceed to present the numerical prescrip- 
tion to construct arbitrary one- and two-qubit gates in 
the next section, it is instructive to first work out some 
important examples whose loop can be constructed an- 
alytically. In particular, we will show that all the gates 
required for the proof of universality may be obtained 
within the present three-state model. 

The first example is the 7r/8-gate, 



U. 



7r/8 



1 

e'""/^ 



(23) 



By inspecting the connection coefficients in Eqs. (|13|-[16D 
we easily find that the loop presented by the sequence 



(02,02) : (0,0)^(7r/2,0)-^(V2,7r/8) 
^(0,^/8)^(0,0). 



yields the desired gate. Note that the loop is in the 
(6*2, 02)-plane and that all the other parameters are fixed 
at zero. Explicitly, we verify that 

U^/s = cxp(^^^Jg^^o)exp (1^2 102=^8) 

X exp -^02le2=V2) '^''P ("i ^e2\^,=o) 
= exp(-^ ^02 162=7^/2) • (25) 
The next example is the Hadamard gate 



H 



V2 U -1 



(26) 



(24) 



Instead of constructing H directly, we will rather use the 
decomposition 

H = e~"/2 exp (i^cTz) exp (i^cFy^j . 

It is easy to verify that the holonomy associated with the 
loop 

(02, 0i) : (0,0)^(V2,0)^(V2,/3) 

^(0,/3)^(0,0) (27) 

is exp(i/3(Ty), while that associated with the loop 

(0i,02,<^i) : (0,0,0)^ (V2, 0,0)^ (V2,7r/2,0) 
~.(V2,7r/2,a)-.(7r/2,0,a) 
^ (0,0, a) -> (0,0,0) (28) 

is exp(za(Tz). Here again, the rest of the parameters are 
fixed at zero. Finally, we construct the phase-shift gate 
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FIG. 1: Objective function landscape in 2D. 



e , which is produced by a sequence of two loops. First 
we construct a gate similar to the 5-shift gate using (cf., 
the 7r/8-shift gate) 

: (0,0)^(7r/2,0)^(7r/2,5) 
->(0,J)->(0,0). (29) 

This loop followed by the similar loop in the {62,4)2)- 
space yields the e*''-gate as 

(01,01,^2,02) : 
(0, 0, 0, 0) ^ (0, 0, 7r/2, 0) ^ (0, 0, 7r/2, S) 
-> (0, 0, 0, S) ^ (0, 0, 0, 0) ^ (7r/2, 0, 0, 0) 

(7r/2, S, 0, 0) ^ (0, S, 0, 0) ^ (0, 0, 0, 0). (30) 

Finally, the controUed-phase gate U (O) = 
exp(i0|ll) (11|) can be written as 

{9^,0 ■■ (0,0)^ (^2,0)-. (V2,e) 

-> (0,8) ^(0,0). (31) 



IV. NUMERICAL METHOD 

Now we adopt a systematic approach to actually con- 
structing arbitrary quantum gates. This is the first time 
that arbitrary one- and two-qubit gates are constructed 
in a three-state model that is in a way the simplest 
possible realization for HQC while still maintaining the 
tensor-product structure necessary for exponential speed- 
up. It has not been shown previously how to construct 
the CNOT, let alone the two-qubit Fourier transform in 
a single loop. Hence, we resort to numerical methods. 
Since it is extremely difficult to see which single loop re- 
sults in a given unitary operator, our approach will be 
that of variational calculus. 

We convert the inverse problem, i.e. which loop cor- 
responds to a given unitary operator, to an optimization 
problem. The problem of finding the unitary operator for 
a given loop is straightforward. Keeping the basepoint 



NODE 


Oi 


O2 


<Pi 


02 


begin 














1 


-5.28 


2.04 


0.18 


-0.40 


2 


-0.44 


1.49 


-0.08 


3.70 


3 


-0.70 


-0.27 


-0.11 


2.59 


end 















TABLE I: Loop of Fig. numerically. 



of the holonomy loop fixed, we let the midpoints vary. 
Owing to the 27r-periodicity, the loops can end either in 
the origin or at any point that is modulo 2tt. 

The space of all possible loops is denoted by V. We 
shall restrict the variational task to the space of polyg- 
onal paths Vfc, where k is the number of vertices in the 
path excluding the basepoint. Naturally, we have Vfc C V 
such that we are not guaranteed to find the best possible 
solution among all the loops, but provided that we use 
a good optimization method, we may expect to find the 
best solution in the limited space Vk ■ Since the dimension 
of the variational space increases with k, one is forced to 
use as low a fc as possible. For instance, for one-qubit 
gates the dimension is 4fc. In the case of two-qubit gates 
the dimension is 9fc. Low k appears to be desirable for 
experimental reasons as well. 

Formally, the optimization problem is to find a 7, such 
that 

f{j)^\\U-U4F (32) 

is minimized over all 7 € Vfc. We naturally hope the min- 
imum value to be zero. Here || • \\f is the so-called Frobe- 
nius trace norm defined by \\A\\p — \/Tt (A^A). We 
could employ the well-known conjugate-gradient method 
to solve the task at hand but this method, or any other 
derivative-based method, is not expected to perform well 
in the present problem due to the complicated structure 
of the objective function. Hence we will use the robust 
polytope algorithm |l4|. 

We have plotted a sample 2D section of the optimiza- 
tion space in Fig. |l|. The axes represent two orthogonal 
directions in the optimization space of a certain two-qubit 
gate. The cc-axis was obtained by interpolating between 
two known minima, whereas the y-axis was chosen ran- 
domly. One can readily verify from the figure that the 
optimization task is indeed extremely hard. 

The calculation of the holonomy requires evaluating 
the ordered product in Eq. (||). The method used in 
the numerical algorithm is to simply write the ordered 
product in a finite-difference approximation by consider- 
ing the connection components as being constant over a 
small difference in the parameters ^7^, i.e. 

« exp {-Adjn)SYn) ■ • • exp {-Ahi)Sjl) . (33) 

Throughout the study we used 200 discretization points 
per edge, i.e., n = 200 x (fc -I- 1). 
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-6 -4-2 ' -6 -4 -2 -6 -4 -2 □ 

01 01 01 

FIG. 2: Loop in parameter space that gives the Hadamard 
gate. 




-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

ei 6*1 6»i 

FIG. 3: Loop in parameter space that yields the gate U — 
exp (ija^) exp (i^ay) expia^. 

V. RESULTS 

First we attempted to find a loop that yields the 
Hadamard gate. Using a random initial configuration, we 
obtained the results that are plotted in Fig. 0. The error 
function /(7) had a value smaller than 10 at the nu- 
merical optimum. The plot represents all the possible pro- 
jections on two perpendicular axes (the horizontal axis 
is always 6i) in the four-dimensional space. Note that 
this optimization was carried out in V3, meaning that 
there are three vertices other than the reference point. 
The results do not take advantage of the 2tt periodic- 
ity. We have also included the data points in Table |. It 
is impressive that such a simple control loop yields the 
gate. Furthermore, this is just one implementation of 
the Hadamard gate. It is possible to find many different 
ones. 

Another example of one-qubit gates is given in 
Fig. I and in Table ||. The gate that we tried 
to implement was now chosen arbitrarily to be C/ = 

exp (iyCTz) exp (i^cTj^) exp ifj^. Again, the error was 
well below 10~* at the optimum. We argue that our 
method is capable of finding any one-qubit gate. These 
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-^1 
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-2.03 


1.31 


0.80 


-1.16 


2 


1.21 


1.18 


-2.35 


0.57 


3 


2.54 


0.66 


-0.49 


0.96 


end 















TABLE II: Loop of Fig. || numerically. 




-2 -1 1 2 3 4 '-2 -1 1 2 3 4 

FIG. 4: Loop in parameter space that gives the CNOT gate. 
Here 7cnot G V3 and the error is below 10"^"^. 



results are not very enlightening as such, but should nev- 
ertheless clearly prove the strength of the technique. 

We also found several implementations for two-qubit 
gates. Figure ^ presents the loop 7cnot G V3 that pro- 
duces the CNOT. We observe, however, that again the 
minimization resulted in an accurate solution. The min- 
imization landscape is just as rough in the case of two 
qubits. Now, of course, the dimension of V3 is 24. 

We also found an implementation of the SWAP gate 
given in Fig. ^. 

Finally, it is interesting to observe that even the two- 
qubit quantum Fourier transform can be performed eas- 
ily. The resulting loop is presented in Fig. It is re- 
markable that such a simple single loop yields a two-qubit 
quantum Fourier transform. We used only three vertices 
but were still able to find an acceptable solution. We ar- 
gue that the error can be made arbitrarily small for any 
two-qubit gate. 
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FIG. 5: Loop in parameter space which realizes the SWAP 
gate. Here the error is below 10~^^. In this case the varia- 
tional space is V5. 

VI. DISCUSSION 

The realization of arbitrary one- and two-qubit gates 
in the context of holonomic quantum computation has 
been demonstrated. By restricting the loops in the con- 
trol manifold within a polygon with k vertices, it be- 
comes possible to cast the realization problem to a finite- 
dimensional variational problem. We have shown explic- 
itly that some useful two-qubit gates are realized by a 
single loop. 

A possible improvement of the present scenario would 
be to minimize the length of the path realizing a given 
gate. This can be carried out by introducing an appro- 
priate penalty or barrier function and the Fubini-Study 
metric in the control manifold <CP^. This optimization 



program is under progress and will be reported elsewhere. 
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FIG. 6: Loop 7Fourior. The error is below 10 
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